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Abstract 

The connection between the complex Sine and Sinh-Gordon equations on the complex plane asso- 
ciated with a Weierstrass type system and the possibility of construction of several classes of mul- 
tivortex solutions is discussed in detail. We perform the Painleve test and analyse the possibility 
of deriving the Backlund transformation from the singularity analysis of the complex Sine-Gordon 
equation. We make use of the analysis using the known relations for the Painleve equations to 
construct explicit formulae in terms of the Umemura polynomials which are r-functions for ratio- 
nal solutions of the third Painleve equation. New classes of multivortex solutions of a Weierstrass 
system are obtained through the use of this proposed procedure. Some physical applications 
are mentioned in the area of the vortex Higgs model when the complex Sine-Gordon equation is 
reduced to coupled Riccati equations. 

PACS subject classifications . Primary 02.30Jr; Secondary 02.30Dk. 



I. INTRODUCTION. 

The complex Sine and Sinh-Gordon equations have been of considerable interest recently in 
many areas of mathematical physics. They originally appeared in the reduction of the 0(4) 
nonlinear sigma modeP'^. They have also appeared in a number of other physical contexts, for 
example, in the study of a massless fermion model with chiral symmetry, and also in the study 
of the motion of a vortex filament in an inviscid incompressible fluid^. The equations have been 
found to be completely integrable, and some work on the construction of multi-soliton solutions 
has been carried out^~^. It has also been shown that the complex Sine-Gordon theory may 
be reformulated in terms of the Wess-Zumino-Wittcn action and interpreted as the integrably 
deformed SU {2) /U{l)-coset model, (as in^ and references therein). These studies are based on 
the complex Sine and Sinh-Gordon theory in the (l-l-l)-dimensional Minkowski space-time. These 
models have energy functionals which are related to that of the Ginsburg-Landau energy 

egl^ J[m'+l{i±mr]d'x. (1.1) 

^Prom the physical point of view, this expression constitutes the basis of the phenomenological 
theory of superfluidity^'^ and has appeared in particle physics as welP. The Ginsburg-Landau 
functional is minimized by the Gross-Pitaevski vortices^. These are topological solitons of the 
form 

V^(x, y) = $„(r)e^"^ $n ^ 1, r ^ oo. 

It is the purpose here to study vortex solutions for the equations of interest discussed in this paper. 

Another important physical application is to the area of superconductivity^, where vortex 
solutions play an essential role. The Lagrangian for the superconducting system usually takes the 
form 

L--If',.+ \V,^\'-V{^). (1.2) 

Here, tp plays the role of the Higgs field, F^i, pertains to the electrodynamic term and V{ijj) is the 
scalar or Higgs potential function which is responsible for mass generation in the system. 

In particular, the classical Sine-Gordon equation has relevance to models which are of interest 
to particle theory. The complete integrability of the Gross-Neveu model, in the large N approx- 
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imation, is quite analogous to the complete integrability of the classical Sine-Gordon equation, 
where 1/N plays the role of the coupling constant^'^. For iV = 1, a large coupling case, the model 
reduces to the massless Thirring model, which is scale invariant with anomalous dimensions. 
For N > 1, the theory exhibits non-trivial renormalization group behavior and mass generation 
through dimensional transmutation. 

Another area of recent importance with regard to applications is to the area of liquid crystals 
and membranes^. Fluid membranes may be idealized as two-dimensional surfaces in an aqueous 
solution with each membrane being made up of a double layer of long molecules. The curved fluid 
membrane may be treated as a bending liquid crystal cell with uniaxial molecular order. Various 
physical properties of interest can be calculated in terms of quantities which are directly related 
to the geometry of the surface. For example, in one model for uniaxial liquid crystals, with the 
normal to the membrane n of the liquid crystal, it has been shown^ that the elastic energy of 
curvature per unit area of the membrane is 

F = J gdA, ^A;(ci + C2 - Cof + kciC2. 

Here, ci, C2 are the two principal curvatures of the surface of the membrane, and cq is called the 
spontaneous curvature of the surface. The quantity, F is referred to as the total bending energy of 
the membrane, the constant Cq is related to the asymmetry of the layers. The positive constant k 
is the bending rigidity and k, which could have either sign, is the elastic modulus of the Gaussian 
curvature. The curvature elastic free energy per unit area of the membrane can also be formulated 
rigorously in terms of two-dimensional differential invariants of the surface^. 

This paper is organized as follows. A generalization of a Weierstrass system for inducing two- 
dimensional surfaces in is presented. It is shown that this system is related to the complex 
Sine and Sinh-Gordon equations. Several properties namely the Lax pair, the Painleve property 
and the Backlund transformation following from the singularity analysis, are investigated for the 
complex Sine-Gordon equation and then extended to the Weierstrass system in Sections 3 and 4. 
Multivortex solutions are constructed and the Painleve structure of the associated radial equations 
is studied. The r functions for the rational class of solutions of the Painleve equation P5 are written 
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in terms of Umemura polynomials, and explicit forms of such solutions are given in Section 5. We 
also present a particular class of solutions of the Weierstrass system via the solution of the complex 
Sine-Gordon equation. Section 6 contains examples of solving the Weierstrass system by means 
of solutions of the complex Sine-Gordon equation. 

II. THE GENERALIZED WEIERSTRASS SYSTEM AND ASSOCIATED COM- 
PLEX SINE AND SINH-GORDON EQUATIONS. 

The Gauss-Codazzi equations describing a two-dimensional surface immersed in a three-dimensional 
sphere which is itself again immersed into a four-dimensional Euclidean space has been studied by 
Darboux^°. He investigated the non-linear Dirac-typc system for four complex valued functions 
ipi and ipi, i — 1,2 satisfying the following system of equations 

d(pi = Q2{^i - ^ ), dip2 = Q2V^2, (^•-^) 

2(p2lp2 

Q,^\lP2?±\lPl\^, Q2 = |¥'2p±|¥'lP, 

and its respective complex conjugate equations. The bar denotes complex conjugate, d — d/dz 
and d = d/dz. The above system can be considered as a variant of the Weierstrass representation^^ 
for surfaces immersed in and will refer to it as such. System (2.1) is a nonlinear first order 
system of eight equations, for which eight of sixteen first order derivatives with respect to z or z 
are known in terms of functions and ipi. System (2.1) admits several conservation laws such as 

d{\ni^2)^d{\ui^2), d{\n^^)^d{\n^2), 9(£i^) = a(^). (2.2) 

As a consequence of conservation laws (2.2), there exist four real- valued functions X^{z, z) which 
are defined by 

— In ■02 + In '02 (i^, X"^ = In (^2 <^-2 + In <^2 c?^, 

(2.3) 

X^ = f In ip2^2 dz + In ■02(^2 dz, X'^ = f /'^^ dz + -J^ dz, 

for any contour F in the complex plane which begins at a fixed zq and ends at variable point z. 
The right hand side of (2.3) does not depend on the choice of the curve F, since the differential of 



equations (2.3) are exact ones. Thus, equations (2.1) and (2.2) allow us to identify the real- valued 

functions X^{z,z), i = I,-- - ,4, as the coordinates of a surface immersed in four- dimensional 
Euclidean space. 

In the present paper, we propose a procedure for constructing explicit multivortex solutions 
of system (2.1). These solutions are obtained through the use of a link between the complex Sine 
or Sinh- Gordon equations on the plane and Weierstrass system (2.1). To our knowledge, this 
connection between these two systems has been observed for the first time here. At this point, we 
want to underline that, for the Weierstrass system (2.1), few explicit solutions have been found up 
to now, and this link with complex Sine-Gordon and complex Sinh-Gordon equations allows us to 
construct new classes of solutions explicitly. We subject system (2.1) to several transformations 
in order to simplify its structure. We start by defining two new complex valued functions 

u^^, v = —. 2.4) 

It is easy to show that if the complex functions ipi and (fi, are solutions of the first order system 
(2.1), then the complex- valued functions u and v defined by (2.4) are solutions of the first order 
system of two equations 

du = -{l±\u\'^)v, dv = --{l±\v\'^)u, (2.5) 

and their respective complex conjugate equations. The elimination of one of the functions u or v 
in system (2.5) leads to the complex Sinh-Gordon (CShG) equation when the sign is positive in 
(2.5), and Sine-Gordon (CSG) equation when the sign is negative in (2.5). Thus we get for both 

cases 

u - 1 

ddu^ —dudu+-uil±\u\^)^Q. (2.6) 

1 ± \uY 4 

As it was shown in^, equation (2.6) was derived in the context of the reduction of the 0(4) nonlinear 
sigma model and as well, the reduction of the self-dual Yang-Mills equations and relativistic 
equations^'^^'^^. 

Note that if v tends to one in CSG equations (2.5), then Q2 vanishes and system (2.1) takes 
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the form 



9^^1=^1(^1 + 77^), ^'^2 = Qli'2. (2.7) 



Conversely, if u tends to one in CSG equations (2.5), then Qi vanishes and system (2.1) becomes 

d<pi = Q2{<Pi - )' 9^-^ = Q2<P2- (2.8) 

2(^2 

These hmits characterize the properties of the solutions of system (2.1). 

Now, we discuss a set of conditions which allow the system (2.1) to become a decoupled system 
of equations. 

Proposition 1. Let the complex functions u and v be solutions of system (2.5). Let the 
functions ■^j and ipi be defined in terms of u and v by 

= eu{l ± |Mn-V2 = ± |^;|2)-V2^ 



(2.9) 



V^2 = e(l ± |Mn-i/2, ip2 = e{l±\v\^)-^''', e = ±l. 
Then the general integrals of system (2.1) are given by 

■01 = uA{z)e^, ifi — vB{z)e^, 

1P2 = A{z)e', (p2 = B{z)e', 
where the complex functions A and B satisfy the following conditions. 



(2.10) 



\A\^ = e-^'+'\l ± \u\Y\ \B\^ = e-^'+'\l±\v\Y^. (2.11) 

Proof: Substituting (2.4) into system (2.1), we obtain an overdetermined system for the 
functions ip2 and (p2 of the following form 



d{u'iP2) = (1 ± \u\''M2\\ui;2 + 7^), di;2 = (1 ± kl')|V'2|'V'2, 

— u 

d{v(p2) = (1 ± \V\^)\(P2\^{V(P2 - — ), dip2 = (1 ± |wp)|<^2|V2. 

2(p2 



(2.12) 



Consider the first equation in the first line of (2.12). Expanding the derivative term d{uip2) and 
using (2.5), this equation reduces to the form 

a^2 = (1± |Mn|V'2|¥2. 



Using (2.9) to replace |V'2p in this result, it reduces to 9'02 = ■02- Making use of (2.9), the 
remaining three equations in (2.12) can be treated in a similar way. Thus the initial system (2.1) 
becomes a linear system of the form 

^■02 = ■02, dip2 = </?2- (2-13) 

These two equations can be easily integrated to give ■02 and (^2 as given in (2.10), then (2.4) can 
be used to obtain ■^i and </7i. The results in (2.10) must be consistent with those in (2.9). If we 
calculate the modulus of ■02 and from (2.10) and equate to the modulus calculated from (2.9), 
the conditions (2.11) arc obtained. In fact, equating ipi, ipi in (2.10) to their corresponding forms 
in (2.9), we must also have that 

A{z)e' = e(l ± |k|')-'/' = A{z)e', B{z)e' = e(l ± \v\Y^/^ = B{z)e\ 

A set of differential constraints which must be satisfied can be obtained by substituting (2.9) into 
(2.13), and we find 

{duu + ^uv{l ± \u\'^)) = ^2(1 ± l-^l^), {dvv - ^uv{l ± \vf)) = :i:2(l ± \v\'^). 

QED 

^Prom the computational point of view, it is more convenient to deal with the CShG or CSG 
equations (2.5) than with the original system (2.1). From every solution of CShG or CSG equations 
(2.5), we can integrate a linear system (2.13), and consequently, a very large class of solutions of 
system (2.1) can be found exphcitly by making use of formulae (2.10) and (2.11). 

Using the connection between the CSG equation (2.6) and Weierstrass system (2.1), we discuss 
in the next section in detail the Painleve analysis of the CSG equation which allows us to extend 
this analysis to the Weierstrass system. 

III. PAINLEVE ANALYSIS OF THE COMPLEX SINE-GORDON EQUATION. 

Integrability of the CSG equation (2.6) is confirmed by tests for the Painleve property. We 
perform the classical test of^^ extended to partial differential equations in^^, assuming a solution 
in the form of a Laurent series about an arbitrary singularity manifold F{z^ ^) = and checking 
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compatibility of the resulting recurrence formulae. Detailed discussion of the meaning and validity 
of this test may be found in^°. The test is carried out for the CSG equation (2.6) and for the system 
(2.5). Both versions of the CSG pass the test. For the system (2.5), we discuss the possibilities 
and limitations of the approach to the Backlund transformation^^. 

Both those forms of the CSG contain non-analytic expressions, therefore z and z should be 
treated as two independent variables and extended to two separate complex planes. Similarly, the 
functions u and u need not be complex conjugates of each other, when their independent variables 
are separately extended. We assume also the same holds for v and v. To avoid the misleading 
complex conjugate symbol, we denote z by u by w and v by s, while symbols of the derivatives 
d and d will be replaced by alphabetic subscripts z and respectively. 

In this notation, the original CSG equation (2.6), is cast into polynomial form. It has the form 
of the system 

(1 — uw)uzt + wUzUt + \u{l — uwY — 0, 

(3.1) 

(1 — uw)wzt + uWzWt + \w{l — uwY — 0. 
Both dependent variables, u and w, begin their series with F inverse. The coefficients at this 
power of F may be written as 

uo = 2{FzFtfl^/Q^, Wo = 2{FzFty/'Qo, (3.2) 

where Qo is an arbitrary function of z and t. 

The indices, at which the determinant of the linear recurrence formulae turns zero, are — 1 
(corresponding to arbitrariness of F), (due to arbitrariness of Qo), 1 and 2. Both compatibihty 
conditions at the last two indices are satisfied. As the 2x2 matrices of the respective linear systems 
have rank one at each of these indices, each of the recurrence formulae yields one arbitrary function 
(say, Qi and (52), to complete the set of four first integrals (F, Qq, Qi, Q2) in the general solution 
of (2.6). To conclude that the CSG equation (2.6) has the Painleve property, we should also 
examine the singularity uw —>■ 1. This task is straightforward when working with system (2.5). 

The system (2.5) is just a form of the CSG equation (2.6), so it should pass the same test. 
We describe it in more detail because of the close connection of (2.5) with the Weierstrass system 
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(2.1). 

For the purpose of the "Painlevc test", equation (2.5) is a 4 x 4 system (we choose the minus 
sign for further analysis) given by 

Uz — {1 — wu)v/2 — 0, Wf — {1 — wu)s/2 — 0, 

(3.3) 

vt + {l- sv)u/2^Q, Sz + {1- sv)w/2^Q. 
We may choose —1 as the initial exponents in the Laurent series of either u and w or v and s. 
We perform the calculations for the former choice, which corresponds to our previous analysis of 
the CSG equation (2.6); the latter choice is equivalent to it by exchange of u with v, w with s 
and z with —t. The above-mentioned singularity uw — > 1 corresponds to the latter choice. The 
resulting initial exponents of v and s are both equal to and their Laurent expansions turn into 
the Taylor ones. Coefficients uq and wq are the same as those in (3.2), while 



^^0 



(F,/F,)VVQo, 50 = (F,/F,)V2g„. (3.4) 



At this point, we would hke to mention the following fact. Usually, the classical "Painleve 
test" is not possible when terms of lowest order in F have nonnegative exponents. Such terms 
may lack the property (necessary for the algorithm of^^ and^^) that differentiation decreases their 
order of magnitude by one. However, here the derivatives of v and s do not contribute to the 
dominant terms in the last two equations of (3.3) and that deficiency is meaningless. 

The remaining terms are derived from the recurrence formulae 



nF^Un + {FjQl)Wn + 2FzFtVn = -{Un-l)z + 
1 \— vn— 1 1 \— vn— 1 v~^"— 's— 1 

"2^0 2^k=l '^kVn-k - 2 Z^fc=0 2^1=1 WkUiVn-k-U 

{FtQl)Un + nFtWn + 2F^FtSn = -{Wn-l)t + |Sn-2 

=0 1^1=1 '^k'^l^n-k-li 



2 

ii\ V^""! lie 1 ^— ^n— 1 ^-vn— A;— 1 ^ 

(n - l)FtVn - {F,/Ql)Sn = -{Vn-l)t 
\n—l I 1 v^"— 1 v^"— 1 



(3.5a) 

(3.56) 
(3.5c) 
(3.5(i) 



+ 2'^0 Z^k=l ^kVfi-k + 2 Z^fe=0 2-^1=1 SkUiVn-k-l, 

-{FtQl)Vn + {n- l)FzSn = -{Sn-l)z 
+ 2^0 2^k=l '^kSn-k + 2 Z^fe=0 2^1=1 VkWiSn-k-l, 

where v and s with negative subscripts are both set equal to zero, and have been admitted for 
compact notation. 



The determinant of the system (3.5) reads 

{F,Ftf{n + l)n{n-l){n-2), (3.6) 

Hence, the indices at which it becomes zero, are —1, 0, 1 and 2 as in the previous case of the CSG 
equation (2.6). Similarly, all the compatibility conditions are satisfied, whence we conclude that 
equation (2.5) has the Painleve property. 

This Painleve integr ability obviously extends to the Weierstrass system (2.1) as the Painleve 
property is invariant under the homographic transformation which converts the equations (2.4) to 
(2.1). 

The method of the Laurent expansion may often be extended to deriving the Backlund trans- 
formation and further an explicit integration scheme^^. The usual approach relies on truncation 
of the Laurent series at the term of order F^, which is assumed to satisfy the original equation. 
The truncation usually is possible through appropriate choice of the arbitrary functions (first inte- 
grals). Some extra assumptions on the coefficients and expansion variable may also be necessary. 
A systematic approach to that problem may be found in 

However, the usual method contains assumptions which are too restrictive for application to 
equations (3.3). Namely, the Laurent series of v and s, which begin with the F° terms, reduce 
to a single term each. This means that v and s would not be transformed at all. Moreover, the 
truncation at F^ implies vanishing of terms proportional to F^. This imposes further constraints 
on these variables: from their recurrence equations (3.5) at n = 1 

-{F,/Ql)si^-{vo)t^O: (3.7a) 

-{FtQl)v, = -{so), = 0. (3.76) 

It follows that vo should be independent of t, while sq should be independent of z. However, these 
coefficients are reciprocals of each other (3.4), whence neither of them may depend on z or t. This, 
together with the truncation of the series, reduces u and v to constants. If we denote 

s — k, V — 1/k, k — constant, (3-8) 
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then the original equations (3.3) reduce to a system of coupled Riccati equations 

u,-{l-wu) /{2k) = 0, wt-{l- wu)k/2 = 0, (3.9) 
which may immediately be linearized by substitution 

{2/k)(M)t, w^2k(\n^)„ (3.10) 

to the Helmholtz equation 

= (l/4)vl/. (3.11) 

Obviously, this linearization also yields an (almost trivial) transformation of u and v, a superpo- 
sition principle, and other properties. 

The method of the Laurent expansion starting from eqs. (3.1) is free of those shortcomings. 
However the condition of vanishing of the higher order coefficients results in cumbersome equations, 
which makes the procedure hardly worthwhile in comparison with other methods of integration 
(see next section). 

IV. ON EQUIVALENCE OF TWO FORMS OF THE COMPLEX SINE-GORDON 
EQUATION. 

Throughout this paper we investigate the CSG equation in the form of equations (2.6). Another 
form of the CSG equation was given in^'''^^ as follows 

Let 

u = sin($/2) exp(iQ;). (4.2) 

The polar coordinates $ and a satisfy a system of equations similar to that given by Lund'' 

sinf$/2) - 1 

^^^-2 iT^dada + - sin ^ = 0, (4.3) 

cos3($/2) 4 

d{tan\^/2)da) + d{tan\^/2)da) = 0. (4.4) 
We formulate the following statement for systems (2.6) and (4.1). 
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Proposition 2. Equation (4.1) is transformed into CSG equation (2.6) through the following 

relations 

1_ 1 

^=-2;, T] — —-Z, 5 = — sin$ exp(— i/?), p = sin $ exp(i/?), (4.5) 
where the phase is given by 

(3^ f {[l + tSin^(^z\z)/2)]d'a{z',z)dz' + [l-ta,n^(^z,z')/2)]d'a(z,z')dz'}. (4.6) 

J zo 

The lower limit of integration, zq, is fixed and depends on the initial conditions. 

Proof. Substitution of (4.5) and (4.6) into (4.1) yields the system (4.3) and (4.4). 

Note that both equations, (2.6) and (4.1) depend on their phases a and f3 respectively, through 
their derivatives only, except for linear dependence on factors exp(iQ;) and exp(i/3). Therefore, any 
change of /3, which leaves its derivatives unchanged, for example, one that arises from deformation 
of the integration contour in (4.6), does not affect equivalence of those equations. Moreover, for 
those u, which satisfy the CSG equation (2.6), the integrand is an exact differential and the path 
of integration does not even affect the value of the phase. 

Note also that the form of the equation determining evolution of the phase (4.4) suggests 
integration in terms of an arbitrary potential iplzjz) 

tan^($/2)9Q; = 9V', tsin^(^/2)da ^ -dijj. (4.7) 

However, the potential is not arbitrary since the compatibility condition dda — dda imposes a 
constraint of the form similar to the original phase equation (4.4) 

9(cot^($/2)9V') + 9(cot^($/2)9V') = 0. (4.8) 

Obviously, if a is a solution of (4.4) for a given $ then ip solves the same equation for $ shifted 
by an odd multiple of tt or subtracted from such a multiple. However, this is not a symmetry 
of the CSG equation (2.6), as the signs in the amplitude equation (4.3) are changed by such a 
transformation. Finally, repetition of the transformation brings us back to the original equation. 
A similar property holds for the CSG equation in the form (4.1), for which the phase equation 
may be written as 

d{cos-^ - d{cos = 0, (4.9) 
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where the change of independent variables has already been performed. 

This symmetry has an additional consequence. The transformations (4.5) and (4.6) from 
equations (2.6) to (4.1) have some features of a Backlund transformation. Namely, the definition 
of (5 by means of the contour integral (4.6) is a solution of a coupled pair of differential equations 



a/3 = [1 + tan2($/2)]9Q;, 9/5 = [1 - ts.n^{^ /2)\da. 



(4.10) 



System (4.10) is overdetermined, as the right hand sides of the above equations must satisfy the 
compatibihty condition 99/3 — Bdp. This condition is indeed satisfied as it proves to be equivalent 
to equation (4.4). Thus we obtain (4.4) in two ways: either through direct substitution of (4.10) 
to (4.9), or from the above compatibility condition. 

The transformation defined by (4.5) and (4.6) may be extended to the inverse scattering 
method. The inverse scattering technique for (4.1) was given in^''. Using (4.5) and (4.6), we 
obtain the Lax pair for (2.6) 



^X = -^[Y^,X], dX=^[Y2 + XY,X], 



(4.11) 



where 



Y 



-i 
i 



Yo^2i 



i cos $ sin ^e*'^ 
sin $6"*^^ —i cos $ 



9(sin$e*^) \ 



cos $ 



(4.12) 



9(sin$e-'^) ^ 
V cos$ 

Further extension to the Weierstrass system (2.1) is also possible, by means of the transformation 
(2.4). The matrices Yi and I2 expressed in terms of u become 



( 



Yo 



( %{\- 2ii2e-2^") yi - (1 - 2ii2e-2i")2e^/^ \ 

\^ — yi - (1 - 2K2e-2i")2e-^/3 _ 2^26-2^") y 

9(^1- (l-2M2e-2««)2e^/3) \ 



V 



9(^1 - (1 - 2ii2e-2i«)2e-^/^) 



1 - 2xi2e-2»" 



1 - 2^26-2^" 





(4.13) 
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^Prom (2.4), complex functions u and v are given in terms of ■01, '02 and (pi, (f2, respectively. 
Taking into account that functions u, v satisfy the same equation (2.6), the resulting Lax pair 
(4.11) for Weierstrass system (2.1) can be described by a system of five two by two matrices Yi, 
Y2 in terms of (V'l, '^2) and {(pi, (p2) respectively, and the constant matrix Y . 
V. MULTIVORTEX SOLUTIONS 

At this point, we would like to derive, through the link between first order system (2.1) and 
equations (2.6), a procedure for constructing multivortex solutions in explicit form. Now, we 
concentrate on certain class of multivortex solutions of equations (2.6) in polar coordinates (r, 6) 
on the plane determined by 

u = An(r)e''*^ n e Z. (5.1) 
Equation (2.6), under the assumption (5.1) is reducible to a second order ODE of the form 

By a homographic transformation of the dependent variable 

A„ = ei+i44, z = r, (5.3) 
1 — w{z) 

where c = —i for the CShG system, and c = 1 for the CSG system, respectively, equation (5.2) 
has the structure of the fifth Painleve (P5) equation 

w"^— —w^ + ^ ^(q;w + -) + + 5.4) 

lw\w — 1) z z w z w — 1 

with the coefficients a and (3 parametrized by a number n e Z and 7, 5 fixed as follows 

a = -p^—, 7 = 0, S=-2. (5.5) 



Such reduction to P5 has recently been performed^^. In general, equation P5 is not integrable in 
terms of known classical transcendental functions. However, for specific values of the parameters, 
solutions of equation (5.4) can be reduced to two types of nontranscendental functions, that is, 

to solutions of a Riccati equation with one arbitrary parameter or to three types of rational 
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solutions of equation P5 According to^^, equation (5.4) with coefficients (5.5) can be written 

in equivalent form as a first order system of ODEs 



dp en 2 ,1 



(5.6) 



'^9 o 2 I A 2 , ^ , 2 

z— = —2z + enq — Az p -\ \- pq , 

dz 2 

where p — w/{l — w). The function q{z) satisfies a Painleve type equation of the form 

// ^ '2 (f' -\- Az"^ q q . 2/ri \ / 2 A 2\2^ 

q = ^5 -T-^q 5 -7-^ h , ^ — 16^2; (2e - n) - (g - Az ) \. 

q^-Az^ q^-Az^ z Az^^q^ - Az^y ^ ' ^ 

The function (f' — Az^ has two roots at g = 22;. Using the transformation 

mz) — , q ^ Zz, 

we get that y{z) is also a solution of equation P5 with parameters 

a^-P^ ^ g \ 7 = 0, (5.7) 

Propositions 3 to 5 are special cases studied by V. Gromak^^ (Chapter 12, section 14) concerning 
the fifth Painleve equations with specific parameters. This analysis is used to construct solutions 
to Weierstrass system (2.1). 

Proposition 3. Let w — w{z) be a solution of the fifth Painleve equation P5 (5.4) with 
parameters given by (5.5), such that the function 

^i{w) = zw'-—w^ + 2zw + — y^0, (5.8) 

does not vanish for any n e Z, then the function 

wi^l- (5.9) 
$i(w) 

is a solution of the fifth Painleve equation (5.4) with parameters given by (5.7). 

Proposition 3 establishes the Auto-Backlund transformation (Auto-BT) for equation P5 when 
7 = 0, (5 = — 2 and a = —j3 are parametrized by n e Z. 
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Now, let us discuss the link between equations P5 with different values of the parameter 5, 
namely, 6^0 and 5 = 0. Note that for 6^0, the solutions of equation (5.4) are expressible 
in terms of Bessel functions, whereas for (5 = 0, they can be expressed in terms of Umemura 
polynomials, which is presented below. 

Proposition 4. Let u{z) 7^ be a solution of equation P5 with parameters given by (5.5). 
Then the function 



where f{z) is defined by 



G? 77/ 1 

/(.) = -ln«(.)--(^.)--^), neZ, 
is a solution of equation P5 with parameters 

\ . f3 = S = 0, 7 = -2- (5-11) 

Based on reference^^ and using the result of Proposition 6, we can find in our case the relation 
between equations P3 with jS ^ 

w" = + -(a^w^ + (3) + ^^w^+—, (5.12) 

w w z w 

and P5 with coefficients given by (5.11). Indeed, the third Painleve equation (5.12) can be written 
as a first order system of ODEs 

zw' = {ea — l)w + e'-)zvP' + zv, e = ±1, 

(5.13) 

zwv' = Pw + Sz + {ea — 2)wv + zv^ . 
^Prom system (5.13), the elimination of w gives 

v" - -^—v 2 + - + 2r2^A^ + ^ + - + -2=0. 5.14 

+ z z^[v^ + d) z^ ^ 6 z^ 

By a homographic transformation of the dependent variable and a change of the independent 
variable 

= Z=V^, (5.15) 
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we obtain from (5.14) equation P5 

where A and B are defined as follows 

A = (5^ + A{-5Y/^(3 - da^ - 45 - 2{-5Y'^ea(3 + 4e5a, 

(5.17) 

B = 5a'^- 2{-6y/ha(3 + 4(-(5)i/2/5 + 45 - /^^ - 4e5Q;. 

Proposition 5. Let y — y{z) be a solution of the fifth Painleve equation (5.16) with param- 
eters given by (5.17), such that the function 

r{z) — w' — {ea — 1) e'jw'^ — 1^0, 

z 

does not vanish. Then the function 

S{r) ^l-2r-\V2^), (5.18) 

is a solution of the third Painleve equation (5.12) with paramters 7 7^ and S — —2. 

The T-functions for the rational class of solutions of the Painleve equation P3 can be constructed^^'^^ 
in terms of the Umemura polynomials T„ = Tn{z, I) which are determined by a sequence of poly- 
nomials in z and defined through the recurrence relation 

with initial conditions Tq = Ti = 1. Based on reference we have 

Proposition 6. For the existence of rational solutions of equation P3 of the form 

W(z) - Tn+l{z,l-l)Tn{z,l) 

^^"^-T„+,(.,/)T„(.,/-l)' ^^-2°^ 

where the Umemura polynomials r„ = Tn{z,l) satisfy recurrence relation (5.19), it is necessary 
and sufficient that the parameters of equation P5 satisfy 

Q/ = 4(n + /), /? = 4(n-/), 7 = -5 = 4. 
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Note that from system (5.13) and the transformation (5.15), there is a connection between 
solutions w of the third Painleve equation (5.12) and the solutions y of the fifth Painleve equation 
(5.16), 

_ zw' - [(4e(n + /)-!)+ 4:ezw]w + 2z , , 

^ ~ zw' - [(4e(n + /) - 1) + 4:ezw]w - 2z' ^' ' 

Substituting the rational solutions (5.20) into formula (5.21) and next replacing the w which 
appears in (5.3) by the function u so obtained, we get multivortex solutions of equations (2.6). 
Consequently, by applying Proposition 2 to the so obtained multivortex solution of (2.6), certain 
classes of solutions of system (2.1) can be found. 

Another class of vortex solutions to CSG equations (2.5) can be provided if we define functions 
u and V in the polar form following 

= >l„(r)e^"^ = A„_l(r)e^("-')^ n e Z. (5.22) 

which transforms the CSG system (2.5) into 

" +-An = {l-Al)An-u 



dr r 

(5.23) 

(^^) - ^ + ^^^n-i = (1 - Al_,)Ar,. 
ar r 

When n — 1, the second equation (5.23) is solved by taking = ±1 and then the first 
equation (5.23) becomes a Riccati equation which can be linearized by a Cole- Hopf transformation 
and solved in terms of Bessel functions. The vortex solution (5.22) takes the form 

u = ^e*^ v = e = ±l, (5.24) 
^o(r) 

where Ji is the Bessel function of the first order, that is, Ji = Io{r), and the prime denotes differ- 
entiation with respect to r. Such a reduction of (5.23) has been recently obtained^^. Consequently, 
from transformation (2.4), we get 

ipi = = e^2. (5.25) 

Substituting (5.25) into Weierstrass system (2.1) and solving the resulting equations, we obtain 

ip2 = F{re-'^), (5.26) 
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where F is an arbitrary function of one variable re and the function satisfies the PDE 

Equation (5.27) has a solution of the form 

i;2^g{v{T)e-''), (5.28) 

where the functions g and v satisfy the following ODEs, 

v' +- - R{r)\^Q, g\+\g\^g^Q, A e C, 
r 

and g denotes the derivative of g with respect to s = v{r)e~''^. These two equations can be 
integrated to give the following expressions 

v(r) = - f tR{t) dr, ^ ^ c. (5.29) 
r Jo 9^ 

^Prom equations (5.25), (5.26) and (5.28), we can summarize the results as follows, 

T'(r) \p^^ r Ap"*" r 

^i = fr^e^'9i^ TR{r)dr), ^ gi- / rR{r)dr). 

Io{r) r Jo r Jo ^^^^^ 

where is a function of one variable, which is restricted by relation (5.29). Note that the solution 
for the function u in (5.24) has the form of a scalar field which has appeared in the study of the 
vortex solutions of superconductivity with asymptotic behavior of the radial part of the solution 
going to zero as r goes to zero and constant as r goes to infinity. Consequently, solutions (5.30) 
of the Weierstrass system (2.1) possess similar asymptotic behavior when functions F and g are 
bounded. 

VI. EXAMPLES OF SOLVING THE WEIERSTRASS SYSTEM VIA THE COM- 
PLEX SINE-GORDON EQUATION 

Now we investigate the possibility of generating new multi-soliton solutions by taking products 
of known solutions of the CSG system (2.6). Thus we can formulate the following: 
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Proposition 7. Suppose it is a solution of equation (2.6) with constant modulus |wp = |cp ^ 

1. Suppose also that a complex function w exists which satisfies |wp = 1, and the differential 
constraint equation 

- wIcP - 1 - - 

u(ddw ^ -j^ _[_ I |2 {dw){dw)) + ^ _^ ^ {dudw + dudw) = 0. (6.1) 

Then the product function U — uw is a, solution of system (2.6). 

Proof: Differentiating the function U = uw, we get the following expressions 

d{uw) — {du)w + u{dw), d{uw) — {du)w + {dw)u, 

dd{uw) — {ddu)w + {du){dw) + {du){dw) + u{ddw). 
Substituting U into equation (2.6), using \u\'^ = |cp and \w\'^ = 1, we obtain that 

{ddu)w + {du){dw) + {du){dw) + u{ddw) 

TiiTi — — mil 

IF j^^((aK)w + u{dw)){{du)w + u{dw)) + ^(1 ± |c|2) (6.2) 

= {ddu)w + {du){dw) + {du){dw) + u{ddw) 

1 ± |cp 1 ± |cp 1 ± |cp 1 ± |cp 4 

Substituting the second derivative ddu from equation (2.6) into equation (6.2) and next collecting 

terms with respect to first derivatives of u and w and simplifying, we obtain the differential 

constraint (6.1). QED 

In the case of the CShG equation, the constant c need not necessarily have modulus different 
from one, since there is no singularity in the (1 + |cp)~^ term in this case. 

At this point, we would like to illustrate Proposition 10 for constructing a solution to system 
(2.6) with an elementary example. The simplest solution of analytic type, a vacuum solution, is 
given by 

u^ce^^--A^)^ (6.3) 

where c and A are complex constants. By substituting the solution u into (2.6), it is easy to show 
that this is a solution provided that the following constraint holds between c and A, 

2e|A| = l±|c|^ e = ±l. 
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Suppose that / is a complex valued function of one complex variable z and define the function w 
as follows, 

/(?)' 

such that f{z) satisfies the constraint (6.1), namely, 

df{z)df{z) + Adf{z)f{z) + Af{z)df{z) = 0. 
Then Proposition 7 implies that the function 

/»' 

is also a solution to system (2.6) and represents a one-soliton solution. 

Let us illustrate our considerations by an example, of the exponential solution (6.3) for A = 
—ia. We show that we can use this to obtain a solution to system (2.1). First of all, since 
eliminating u from the pair of equations in (2.5) results in an equation which is identical to (2.6) 
but with u replaced hj v, we can assign the solution obtained from the second order equation to 
either the u or the v variable. Let us take u — ce*"^""'"^-* . The second function v is obtained from 
(2.5), that is 

2(9?/ 

u = ce^"("+^), V = = ice'^^'+'l (6.4) 

1 ± \uf 

These satisfy (2.5) provided that 2a = 1 ± jc^. Prom (2.4) we can write ■^i = ^■02 and cpi — v(p2, 
and calculate the quantities (5i,2 from (2.1) 

Q, = |^2r(l ± |c|'), Q2 = 1^2^(1 ± |c|'). (6.5) 

Now, we will show that we can find a particular class of solutions which satisfy (2.1). Pirst from 
equations (2.1) with (6.5), the equation 

a^2 = |^2r(l±|c|') 

is satisfied by a function of the form ^2 = e'^^^^^^ provided that b = —2a to be consistent with 
the condition 2a = 1 ± |cp. Similarly, a function of the form ip2 — e^"^^"^-' satisfies the equations 
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dip2 — The equations for ■^i and </?i in (2.1) can also be integrated, and the set of functions 

below 

(6.6) 

constitutes a solution to system (2.1). Substituting (6.6) into (2.3) and integrating, we obtain the 
parametric form of a surface. One obtains that 

= -a{z - zf + 2a\z\^ = ay^ + 2ar^, X'^ = a{z - zf - 2a\z\^ = -ay'^ - 2ar^ 

X^ = 0, X^ ^ i\c\\z - z) ^ -\c\^y. (6.7) 

where we set 2; — z = iy and \z\'^ — r^. Treating y and r > as parameters, one can plot X^, X"^ 
and X"^ to give the surface given in Figure 1, which has the form of a parabolic cylinder. Moreover, 
from (6.12), we can calculate the components of the induced metric 

4 

9.Z = = 8a\z^ - 4\z\^ + 4f) - \c\^ = fe, 

i=l 

4 

9z-z = J]^:^.- = 8a\2z + z){z- 2z) + \c\\ 

A procedure for obtaining solutions to (2.1) as a result of using Proposition 10 can be developed 
from the above example. The new product solution can be called either u or v. One substitutes 
this in the corresponding equation in (2.5) to obtain the remaining unknown solution v or u. 
Using these results in (2.4) to eliminate (■02, (^2), one tries to integrate the nonlinear system (2.1) 
to obtain the required complex functions ■^i and </?i, which are solutions of the Weierstrass system 
(2.1). 

IV. FINAL REMARKS. 

Equations (2.1) and (2.5) under investigation here have long been of interest in field theory. 
In particular, there has been the extensive use of applying soliton solutions to construct models 
of extended particles The Sine-Gordon equation, it seems is the only Lorentz- invariant, 

nonlinear equation whose initial value problem has been solved ^. This equation also describes a 
completely integrable Hamiltonian system. It would certainly be of great interest to find other 
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Lorentz invariant integrable systems. Of more recent interest is the study of vortex tubes The 
motion of vortex tubes in an inviscid incompressible fluid is described by the Biot-Savard law. 
The recently proposed localized induction equation is the simplest model to capture the leading 
order behavior of the three-dimensional self- induced motion of a vortex filament. This type of 
equation is in fact related to the cubic nonlinear Schrodinger equation for a complex variable, and 
implies that the localized induction equation is completely integrable. Note that from solutions 
(3.36), when the functions F and g are real polynomials in a single variable, the vortex structure 
of the solutions is preserved at the level of the functions ■0a and ^pp. These functions are applied 
to generate surfaces. Consequently, by plotting such results in two or three dimensions, these 
surfaces could model a vortex filament in such a fluid 

We have presented a new approach to the study of the Weierstrass system (2.1) in connection 
with CShG and CSG equations (2.6). It proved to be particularly effective in constructing multi- 
vortex solutions of (2.1) in terms of r-f unctions based on rational solutions of the third and fifth 
Painleve equations. It is worth noting that the approach to the Weierstrass system (2.1) proposed 
here can be applied with some necessary modifications, to more general cases of Weierstrass type 
systems describing more diverse surfaces immersed in multi-dimensional Minkowski and pseudo- 
Riemannian spaces. The task of obtaining new types of minimal surfaces described by system 
(2.1) will be undertaken in our future work. 
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Figure Captions. 

Figure 1. The surface which corresponds to the solution (6.6) of the CShG equation and 

coordinates (6.7) for the choice of constants a = 1 and c = 1. 
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